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ABSTRACT 

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a deep z'-band imaging 
survey covering the Spitzer SWIRE Legacy fields designed to create the first large homogeneously- 
selected sample of massive clusters at z > 1 using an infrared adaptation of the cluster red-sequence 
method. We present an overview of the northern component of the survey which has been observed 
with CFHT/MegaCam and covers 28.3 deg^. The southern component of the survey was observed with 
CTIO/MOSAICII, covers 13.6 deg^, and is summarized in a companion paper by Wilson et al. (2008). 
We also present spectroscopic confirmation of two rich cluster candidates at z ~ 1.2. Based on Nod- 
and-Shuffle spectroscopy from GMOS-N on Gemini there are 17 and 28 confirmed cluster members in 
SpARCS J163435+402151 and SpARCS J163852+403843 which have spectroscopic redshifts of 1.1798 
and 1.1963, respectively. The clusters have velocity dispersions of 490 ± 140 km/s and 650 ± 160 
km/s, respectively which imply masses (M200) of (1.0 ± 0.9) x lO" Mg and (2.4 ± 1.8) x 10^'' Mq. 
Confirmation of these candidates as bona fide massive clusters demonstrates that two-filter imaging 
is an effective, yet observationally efficient, method for selecting clusters at z > 1. 

Subject headings: infrared: galaxies 



1. INTRODUCTION 

In the nearby universe there are numerous lines of ev- 
idence suggesting that environmental processes could be 
the dominant force driving the evolution of the galaxy 
population. Properties such as star formation rate (SFR, 
e.g., Lewis et al. 2002; Gomez et al. 2003; Kauffmann 
et al. 2004), morphology (e.g.. Dressier 1980; Goto et al. 
2003; Park et al. 2007), stellar mass (e.g., Kauffmann 
et al. 2004), color (e.g., Hogg et al. 2003; Balogh et al. 
2004; Blanton et al. 2005), and luminosity (e.g., Croton 
et al. 2005; Park et al. 2007) are all strongly correlated 
with local galaxy density. Although there is still debate 
about which, if any, of these relations are "fundamental" 
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(e.g., Hogg et al. 2004; Park et al. 2007), it is clear 
that the mean properties of galaxies we measure depend 
strongly on the type of environment they occupy. 

An obvious first step toward a better understanding of 
environmental processes is to study how they evolve with 
redshift. At higher redshift, the overall population of 
galaxies is younger and have been living within their local 
environment for less time. The environmental processes 
that are most effective and have the shortest timescales 
should be most apparent when comparing galaxies at dif- 
ferent densities in the high redshift universe. The data at 
higher redshift are still somewhat sparse compared to the 
nearby universe but it is beginning to emerge that prop- 
erties such as the SFR (e.g., Elbaz et al. 2007; Cooper 
et al. 2008; Poggianti et al. 2008), color (e.g.. Cooper 
et al. 2007) and morphology (e.g.. Dressier et al. 1997; 
Postman et al. 2005; Smith et al. 2005, Capak et al. 
2007) are still correlated with local density, albeit differ- 
ently from the nearby universe. 

Of particular interest for understanding environmen- 
tal processes are the cores of rich galaxy clusters. These 
are the most extreme density environments at all red- 
shifts, and if environment is truly an important force in 
galaxy evolution a comparison of the properties of galax- 
ies that live in this environment to those that live in the 
field should provide the largest contrasts. Despite their 
potential value for such studies, and the abundance of 
resources directed at finding distant clusters, there are 
still relatively few confirmed rich clusters at z > 1. 

The major challenge for cluster surveys targeting the 
z > 1 range is the need to be simultaneously deep enough 
to detect either the galaxies or hot X-ray gas in clusters 
and yet wide enough to be able to cover a large area be- 
cause of the rarity of rich clusters at z > 1. The require- 
ment of both depth and area has pushed X-ray detection 
of clusters with current telescopes to the limit. 
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The largest area targeted X-ray eluster surveys are the 
XMM-LSS (Valtchanov et al. 2004; Andreon et al. 2005; 
Pierre et al. 2006) and the XMM-COSMOS (Finoguenov 
et al. 2007), and while these have been successful at dis- 
covering z > 1 clusters (e.g., Bremer et al. 2006), they 
cover areas of only 9 and 2 deg^, respectively and are 
therefore limited to fairly low mass systems on average. 
Indeed, X-ray detection of clusters at > 1 is so chal- 
lenging that currently the most promising surveys are 
those searching for clusters serendipitously in the entire 
XMM-Newton archive (e.g., Romer et al. 2001; Mullis 
et al. 2005; Stanford et al. 2006; Lamer et al. 2008). 

Complementary to X-ray detection is optical detec- 
tion of clusters using overdensities of galaxies selected 
using the red-sequence (e.g., Gladders fc Yee 2000, 2005; 
Gilbank et al. 2004; Muzzin et al. 2008) or photomet- 
ric redshifts (e.g., Stanford et al. 2005; van Breuklen et 
al. 2007; Eisenhardt et al. 2008). Recently, it has be- 
come clear that the key to discovering clusters above z > 
1 with these techniques is the incorporation of infrared 
(IR) data which probes the peak of the stellar emission 
for galaxies at z > 1. 

Although IR surveys have thus far been confined to 
modest areas (ranging from 0.5 - 8.5 deg^) they have 
been extremely successful at detecting z > 1 clusters 
(e.g., Stanford et al. 2005; Brodwin et al. 2006; van 
Breuklen et al. 2007; Zatloukal et al. 2007; Krick et 
al. 2008; Eisenhardt et al. 2008; Muzzin et al. 2008). 
The IR cluster community is now regularly discovering 
clusters at 2: > 1 and shortly IR-detected clusters should 
outnumber their X-ray counterparts. 

Currently the largest area IR survey still deep enough 
to detect clusters at 2; > 1 is the Spitzer Wide-Area 
Infrared Extragalactic Survey (SWIRE, Lonsdale et al. 
2003; Surace et al. 2005). SWIRE covers ~ 50 deg^ in 
the Spitzer bandpasses and is slightly deeper, and nearly 
a factor of six larger than the next largest IR cluster 
survey, the IRAC Shallow Survey Cluster Search (ISCS, 
Eisenhardt et al. 2008). 

In Muzzin et al. (2008) we demonstrated the poten- 
tial of using the red-sequence method with Spitzer data 
to detect distant clusters using data from the 3.8 deg^ 
Spitzer First Look Survey (FLS, Lacy et al. 2004). In 
2006 we began observations for the Spitzer Adaptation 
of the Red-sequence Cluster Survey (SpARCS), a deep z'- 
band imaging survey of the SWIRE fields. SpARCS aims 
to discover the first large, yet homogeneously-selected 
sample of rich clusters at z >1 using the red-sequence 
method. SpARCS is similar to the RCS surveys (Glad- 
ders & Yee 2005; Yee et al. 2007) which target clusters 
to z ~ 1 using an R - z' color except that we use a z' 
- 3.6/xm color, which spans the 4000 A break at 2; > 1. 
With a total area effective area of 41.9 deg^ SpARCS is 
currently the only z > 1 cluster survey that can discover 
a significant number of rare rich clusters. These clusters 
will be extremely valuable for quantifying the evolution 
of galaxy properties in the densest environments at high 
redshift. 

This paper is organized as follows. In § 2 we provide a 
brief overview of the northern component of the SpARCS 
survey (the southern component is summarized in Wil- 
son et al. 2008). In § 3 we discuss the selection of cluster 

candidates that were chosen for foUowup spectroscopy, 
and in § 4 we present spectroscopic confirmation of two 



z > 1 clusters from early SpARCS data. In § 5 we present 
the dynamical analysis of the clusters followed by a dis- 
cussion of the cluster properties in § 6. We conclude with 
a summary in § 7. 

Throughout this paper we assume an n„i = 0.3, Ha = 
0.7, Hq = 70 km s^^ Mpc^^ cosmology. All magnitudes 
are on the Vega system unless indicated otherwise. 

2. THE SPARCS-NORTH SURVEY 

The SWIRE survey is located in six fields and con- 
tains ~ 50 deg^ of imaging in the four IRAC bandpasses 
(3.6/im, 4.5/im, 5.8/zm, and 8.0/im) and the three MIPS 
bandpasses (24/im, 70/im, and 160/im). Three of the 
fields are located in the northern hemisphere (ELAIS- 
Nl, ELAIS-N2, and the Lockman Hole), two of the 
fields are located in the southern hemisphere (ELAIS- 
Sl, and the Chandra-S), and one of the fields is equa- 
torial, the XMM-LSS field. A thorough discussion of 
the data reduction, photometry, cluster finding, and the 
SpARCS catalogue for all fields will be presented in a 
future paper by A. Muzzin et al. (2009, in prepara- 
tion) . Here we present a brief summary of the z'-band ob- 
servations of the ELAIS-Nl, ELAIS-N2, Lockman Hole, 
and XMM-LSSi'^ fields obtained with CFHT/MegaCam; 
hereafter the SpARCS-North Survey. Observations of 
the ELAIS-Sl and Chandra-S fields were obtained with 
the CTIO/MOSAICII and are outlined in the companion 
paper by Wilson et al. (2008). 

The IRAC imaging of the ELAIS-Nl, ELAIS-N2, Lock- 
man, and XMM-LSS fields covers areas of 9.8, 4.5, 11.6, 
and 9.4 deg^, respectively. In Figure 1 we plot the IRAC 
3.6^m mosaics for these fields. The superposed white 
squares represent the locations of the CFHT /MegaCam 
pointings. The pointings were designed to maximize the 
overlap with the IRAC data, but to minimize the over- 
all number of pointings by omitting regions that have 
little overlap with the IRAC data. There are a total of 
12, 5, 15, and 13 MegaCam pointings in the ELAIS-Nl, 
ELAIS-N2, Lockman Hole and XMM-LSS fields, respec- 
tively. 

We obtained observations in the z'-band with 6000s 
of integration time for each pointing in the ELAIS-Nl, 
ELAIS-N2 and Lockman Hole fields in queue mode using 
CFHT/MegaCam which is composed of 36 4096 x 2048 
pixel CCDs, and has a field of view (FOV) of ~ 1 deg^. 
Omitting the large chip gap areas and regions contami- 
nated by bright stars, the total overlap region with both 
z' and IRAC data is 28.3 deg^. 

Photometry was performed on both the z' and 
IRAC mosaics using the SExtractor photometry package 
(Bertin & Arnouts 1996). Colors were determined us- 
ing 3 IRAC pixel (3.66") diameter apertures. The IRAC 
data was corrected for flux lost outside this aperture due 
to the wings of the PSF using aperture corrections mea- 
sured by Lacy et al. (2005). Total magnitudes for the 
IRAC photometry were computed using the method out- 
lined in Lacy et al. (2005) and Muzzin et al. (2008). The 
5(7 depth of the z' data varies depending on the seeing 
and the sky background; however, the mean 5a depth for 
extended sources is z' ~ 23.7 Vega (24.2 AB). 

3. CLUSTER SELECTION 

1* The XMM-LSS data was obtained as part of the CFHT 
Legacy Survey 
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Fig. 1. — The 3.6^tm mosaics for the four SWIRE fields observable from the northern hemisphere. The location of the SpARCS z'-band 
CFHT/MegaCam pointings are overplotted as white boxes. Each MegaCam pointing covers ~ 1 deg'^. The observations of the XMM-LSS 
field were obtained as part of the CFHTLS-wide survey. Excluding areas masked by bright stars and missed by MegaCam chip gaps there 
are 28.3 deg^ with both z' and 3.6/im observations in the northern fields that can be used for cluster finding. 



Clusters are found in the data using the cluster red- 
sequence algorithm developed by Gladders & Yee (2000; 
2005). Muzzin et al. (2008) used a sHghtly modified ver- 
sion of the algorithm to detect clusters at < z < 1.3 in 
the FLS using an R - 3.6/xm color. We use the Muzzin et 
al. (2008) code for the SpARCS data. The change from a 
R - 3.6/im color to a z' - 3.6/im is optimum for targeting 
clusters at z > 1, where the z' - 3.6/im color spans the 
4000A break. Other than the change of using a different 
optical band, the SpARCS algorithm is identical to that 
presented in Muzzin et al. (2008) and we refer to that 
paper for further details of the cluster finding technique. 

After the first semester of z' observations were 
complete there were ^ 14 deg'^ of data with both 
z' and 3.6/im data. From this area we selected 



two rich cluster candidates, both from the ELAIS- 
N2 field with red-sequence photometric redshifts-'^^ 
of z > 1.2 for spectroscopic foUowup. These 
two cluster candidates, SpARCS J163435+402151 
(R.A.: 16:34:35.0, Decl:-|-40:21:51.0), and SpARCS 
J163852-I-403843 (R.A.:16:38:52.0, Decl:+40:38:43.0), 
have richnesses, parameterized by Bgc,_R, of 1053 ± 278 
Mpc^-* and 988 ± 270 Mpc^-^, respectively. For a dis- 
cussion of Bgc and Bgc _r, as cluster richness estimates see 
Yee & Lopez-Cruz (1999), and Gladders & Yee (2005). 
Based on the empirical calibration of Bgc vs. M200 de- 
termined by Muzzin et al. (2007) in the K-band for the 
CNOCl clusters at z 0.3, this implies M200 = 5.7 x 10^'' 



Based on a. Zf 
model 



2.8 single-burst Bruzual & Chariot (2003) 



4 



S 40 

S 20 



: SpARCS_J16345+4021_ml_3 

^ z = UT9 

-. ?: (Vega) =^ 21.54 

— :BL3560 


[oni 

H9:CN HE K H 




G Hj' ^ 


~ SpARCS_J16J4.5 + 4021_ml_21 
'- z' IVega) = 22.14 

— BL3560 


HS 04 K H 




J Uy 


~ SpARCS_J1634.5 + 4021_ml_22 

~- /■ (Vega) = 21.92 

BL3560 


[OH] J. n 

H0:cN h; K H 


Ha 


O tfy -_ 


~ SpARCS.J1634,5 + 4021_ml_23 

- I = 1.180 

'- I'(V^) = 21.45 

^ BL356I} 


ion] 

He:cN HE :k H 




A i 

(! :H;- - 


L SpARCS_J16J4.S+4021_ml_37 

: z = 1.180 

--_ (Vega) = Z1.H2 

— :B13560 


[Oil] /L 
He:cN :he :k h 


II 




r SpARCS_J1634.5 + 4021_iii3_2 

i- z = 1.175 

E- z' (Vega) = 20.51 

L BL3560 


oil] 

HB CN U5 K H 




3 Hi' J 


— SpARCS.J1634.5 + 4021_iii3_69 

: z = 1.177 

-_ z' (Vega) = aj.95 


[OII| 






- :BL3560 


He:cN h; :k :h 






Sl^Spectnim 





SpARCS J163a.a+4039_ml_5 
= 1.194 
(Vega) = 21.52 



SOOO S500 9000 

Wavelength (Angstroms) 



Fig. 2. — Spectra for a subsample of seven galaxies in the clus- 
ter SpARCS J163435-I-402151. The spectra have been smoothed 
with a 7-pixel boxcar so the sampling matches the instrumental 
resolution. 

M0 and 5.1 x lO" Mq for SpARCS J163435+402151 
and SpARCS J163852+403843, respectively. Although 

Muzzin ct al. (2007) found there is a fairly large scatter 
in the Bgc vs. M200 relation, the high richnesses im- 
ply that these candidates are likely to be massive, high- 
redshift systems. 

4. SPECTROSCOPIC DATA 

Multislit nod-and-shuffle (N&S) spectroscopy of 
galaxies in SpARCS J163435+402151, and SpARCS 
J163852-I-403843 were obtained using GMOS-N on Gem- 
ini as part of the program GN-2007A-Q-17. We used 
the R150 grating blazed at 7170A with 1" width slits. 
This provided a resolving power of R = 631 which corre- 
sponds to a resolution of ~ llA, or ~ 250 km s~^ at the 
estimated redshift of the clusters. For all observations 
we used 3" long microslits, corresponding to roughly 
four seeing-disks, allowing a two seeing-disk spacing be- 
tween the nod positions. We observed three masks for 
SpARCS J163435+402151 and four masks for SpARCS 
J163852-I-403843. One mask for each of the clusters was 
observed in "micro-shuffle" mode, but the majority were 
observed in "band-shuffle" mode. All masks were ob- 
served using the RG615 filter which blocks light blueward 
of 6I50A so that multiple tiers of slits could be used. 

Unlike micro-shuffle where the shuffled charge is stored 
directly adjacent to the slit location, band-shuffle shuffles 
the charge to the top and bottom third of the chips for 
storage. While technically it is the least-efficient N&S 
mode in terms of usable area for observations (only the 
central 1.7' of the total 5' FOV can be used) it is ex- 
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Fig. 3. — Same as Figure 2 but for galaxies in the cluster SpARCS 
J163852-I-403843. 

trcmcly efficient for observations of high-redshift clusters 
because it allows the microslits to be packed directly be- 
side each other in the cluster core with no requirement for 
additional space for storing the shuffled charge. In band- 
shuffle mode we were typically able to locate between 
20-26 slits, including three alignment stars, per mask in 
the central 1.7' around the cluster. At 2; ~ 1.2 the 1.7' 
FOV corresponds to a diameter of 850 kpc, roughly the 
projected size of a massive cluster. 

Slits were placed on galaxies with priorities in the fol- 
lowing order: Priority 1, galaxies with colors ± 0.6 mag 
from the red-sequence and 3.6/im < 16.9. Priority 2, 
galaxies with colors ± 0.6 mag from the red-sequence, 
3.6/xm > 16.9 and z' < 23.5. Priority 3, galaxies with 
colors > 0.6 bluer than the red-sequence, but < 1.0 mag 
bluer and 3.6/im > 16.9 and z' < 23.5. Priority 4, same 
as priority 3 but for galaxies with colors bluer than the 
red-sequence by 1.0-1.4 mag. Priority 5: all galaxies 
with 23.5 < z' < 24.5. Roughly speaking, Priorities 1 
through 4 can be described as bright red-sequence, faint 
red-sequence, blue cloud, and extreme blue cloud galax- 
ies, respectively. 

For each mask we obtained a total of 3 hrs of inte- 
gration time by combining six exposures with 30 mins 
of integration time. The six frames were obtained using 
15 nod cycles of 60s integration time per cycle. Each of 
the 6 exposures was offset by a few arcseconds using the 
on-chip dithering option. 

4.1. Data Reduction 

Data were reduced using the CMOS IRAF package. 
We subtracted a bias and N&S dark from each frame. 



The N&S darks arc taken using the same exposure times 
and using the same charge shuffling routine as the sci- 
ence observations, but with the shutter closed. Regions 
with poor charge transfer efficiency cause electrons to be- 
come trapped during the repeated charge shuffling used 
in the observations. Such charge traps can be identified 
and corrected using dark frames taken with the same 
N&S settings. Images were registered using bright sky 
lines and sky subtracted using the complementary stor- 
age area using the "gnsskysub" task. Final mosaics are 
made by coadding the sky subtracted images. 

One dimensional spectra were extracted using the 
iGDDS software (Abraham et al. 2004). Wavelength 
calibration for each extracted spectrum was performed 
using bright sky lines from the unsubtr acted image, also 
with the iGDDS software. Wavelength solutions typi- 
cally have an rms < O.SA. We determined a relative flux 
calibration curve using a long slit observation of the stan- 
dard star EG131. 

Redshifts were determined interactively for each spec- 
trum by comparing with the templates available in 
iGDDS. Most of the redshifts were identified using the 
early-, intermediate-, and late-type composite spectra 
from the Gemini Deep Deep Survey (Abraham et al. 
2004). The final redshifts was determined using the av- 
erage redshift from all absorption and emission lines that 
were detected. The vast majority of redshifts were de- 
termined by identifying the [Oil] 3727A-doublct emission 
line (which is not resolved at our resolution), or the Cal- 
cium II H-l-K absorption lines. Many of the spectra also 
show the Balmer series lines. We list the spectroscopic 
members of SpARCS J163435+402151 and SpARCS 
J163852-h403843 in Tables 1 and 2, and the spectroscop- 
ically confirmed foreground/background galaxies in Ta- 
bles 3 and 4. We also plot examples of some cluster 
galaxy spectra in Figures 2 and 3. R, z', and 3.6/xm 
color composites of the two clusters are shown in Figures 
4 and 5. The white squares denote the spectroscopically 
confirmed cluster members and the green squares denote 
the spectroscopically confirmed foreground/background 
galaxies. 



5. CLUSTER VELOCITY DISPERSIONS 

For both clusters there are a sufficient number of red- 
shifts to determine a velocity dispersion ((t„) and there- 
fore a dynamical mass. Our cr^'s are determined us- 
ing the method detailed in Blindert (2006). Briefiy, we 
make a rejection of near-field non-cluster members using 
a modified version of the Fadda et al. (1996) shifting- 
gap procedure. This method uses both the position and 
velocity of galaxies to reject interlopers. In Figure 6 we 
plot the relative velocities of the cluster galaxies as a 
function of projected radius. Two galaxies in SpARCS 
J163435H-402151 are rejected as near-field interlopers 
and three are rejected in SpARCS .J1G3852-F403843. The 
rejected galaxies are plotted as crosses in Figure 6 and 
are not used in computing the mean redshift of the clus- 
ter or ay. Once outliers are rejected the redshift of 
the clusters is determined using the remaining galax- 
ies. The spectroscopic redshift of the clusters is 1.1798 
and 1.1963 for SpARCS J163435-h40215 and SpARCS 
J163852-F403843, respectively. 



TABLE 1 

Spectroscopic Cluster Members in SpARCS 
J163435-I-402151 



ID R.A. Decl. z' Zsp. 

J2000 (Deg.) J2000 (Deg.) Mag Vega 



Mask 1 


3 


248.6589 


40.35303 


21.54 


1.179 


5 


248.6497 


40.36148 


21.58 


1.181 


6 


248.6467 


40.36418 


20.95 


1.166 


21 


248.6513 


40.35238 


22.14 


1.174 


22 


248.6605 


40.35474 


21.91 


1.178 


25 


248.6510 


40.35855 


21.44 


1.180 


28 


248.6605 


40.36132 


22.19 


1.182 


34 


248.6527 


40.36477 


22.78 


1.182 


37 


248.6496 


40.36696 


21.82 


1.180 


41 


248.6459 


40.36841 


21.97 


1.185 


45 


248.6440 


40.37413 


22.14 


1.187 


3027 


248.6546 


40.35461 


21.59 


1.181 


Mask 2 


20 


248.5991 


40.34870 


23.10 


1.170 


27 


248.6548 


40.35104 


22.23 


1.164 


36 


248.6171 


40.36649 


22.36 


1.184 


57 


248.6700 


40.39159 


22.80 


1.176 


Mask 3 


2 


248.6584 


40.34934 


20.50 


1.175 


38 


248.6474 


40.36152 


21.67 


1.178 


69 


248.6669 


40.39132 


21.89 


1.178 



TABLE 2 

Spectroscopic Cluster Members in SpARCS 
J163852-I-403843 



ID R.A. Dccl. z' z, 

J2000 (Deg.) J2000 (Deg.) Mag Vega 



Mask 1 


5 


249.7132 


40.63952 


21.51 


1.194 


7 


249.7144 


40.64159 


21.68 


1.194 


8 


249.7152 


40.64527 


21.41 


1.195 


1007 


249.6974 


40.63383 


22.50 


1.199 


1010 


249.7151 


40.63804 


22.48 


1.202 


1016 


249.7006 


40.64276 


22.20 


1.200 


1025 


249.7137 


40.64952 


22.46 


1.200 


1028 


249.7104 


40.65532 


21.96 


1.190 


1031 


249.7029 


40.65968 


22.85 


1.176 


3026 


249.6668 


40.64420 


22.73 


1.195 


Mask 2 


1018 


249.7032 


40.64533 


22.67 


1.186 


1020 


249.7099 


40.64631 


22.39 


1.202 


1024 


249.7135 


40.64828 


22.24 


1.198 


1026 


249.7126 


40.65311 


22.26 


1.200 


1030 


249.7581 
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Fig. 4. — Left: Rz'SSfim color composite of the cluster SpARCS J163435+402151 at z = 1.1798. The R and z' images have been 
convolved to match the 3.6/xm PSF. The FOV of the image is ~ 3.5' across. Right: Same as left panel but with spectroscopically confirmed 
cluster members marked as white squares and spectroscopically confirmed foreground/background galaxies marked as green squares. 
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Fig. 5.— As Figure 4, but for the cluster SpARCS J163852+403843 at z = 1.1963. The FOV of the images is ~ 4.5' across. 
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timator for systems with > 15 members, and the gap- 
per estimator for systems with < 15 members. As dis- 
cussed in those papers and Blindert (2006) these estima- 
tors are more robust than standard deviations as they 
are less sensitive to outliers, which may still persist even 
after the initial shifting-gap rejection. Using the "ro- 
bust" estimator, SpARCS J163435-h402151 and SpARCS 
J163852+403843 have cr„ = 490 ± 140 km s"! and 650 
± 160 km s~^, respectively, where the errors have been 
determined using Jackknife resampling of the data. 

We estimate the dynamical mass using M200: the mass 
contained within r2ooj the radius at which the mean in- 
terior density is 200 times the critical density [pc)- We 
use the equation, 

M200 = ^Trrfoo • 200pc, (1) 

with the dynamical estimate of r2oo from Carlberg et al. 

(1997), 

where H(z) is the Hubble constant at the redshift of 
the cluster. From these relation we derive r2oo = 
0.62 ± 0.18 Mpc, and 0.82 ± 0.20 Mpc for SpARCS 
J163435-I-402151 and SpARCS J163852+403843, respec- 
tively. From Equation 2, these imply M200 = (1-0 ± 
0.9) X 10" Mq and (2.4 ± 1.8) x lO" Mq for SpARCS 
J163435-F402151 and SpARCS J163852-h403843, respec- 
tively. 

6. DISCUSSION 

6.1. Red-sequence Photometric Redshifts 

In Figure 7 we plot the z' - 3.6/um vs. 3.6/um color mag- 
nitude relation for galaxies at projected radii (R) < 550 
kpc in the fields of both clusters. Spectroscopically con- 
firmed members and confirmed foreground/background 
galaxies are plotted as red and blue diamonds, respec- 
tively. The dotted line in both panels is the best fit 
line with the slope fixed at zero, to the spectroscop- 
ically confirmed members. These lines indicate that 
the red-sequence galaxies have z' - 3.6/im colors of 4.77 
and 4.82 for SpARCS J163435-h402151 and SpARCS 
J 163852-1-403843, respectively. Using a solar metallic- 
ity, Bruzual & Chariot (2003) simple stellar population 
(SSP) with a. zf = 4.0 these colors imply photometric 
redshifts of 1.19 and 1.20, in excellent agreement with 
the spectroscopic redshifts. At z ^ 1.2, the red-sequence 
photometric redshifts do not depend strongly on the cho- 
sen z/. If we instead use a z/ = 2.8 SSP, the red-sequence 
color would predict redshifts of 1.21 and 1.24, and for a 
Zf = 10.0 SSP it would predict redshifts of 1.13 and 
1.15. However, the color differences between all these 
models at fixed redshift are small (< 0.1 mag), and so it 
is not possible to distinguish between different formation 
epochs without more data. Still, the close agreement 
between the red-sequence photometric redshift derived 
using a reasonable zj and the spectroscopic redshift is 
encouraging for the use of red-sequence photometric red- 
shifts for clusters at z > 1. 

6.2. Mass vs. Richness 

Both clusters have lower masses than predicted by 
their richness by factors of ~ 6 and 2 for SpARCS 



J163435+402151 and SpARCS J163852-h403843, respec- 
tively, although due to the large error bars the differ- 
ences are only significant at ~ 1 and 2a, respectively. 
Whether this represents a redshift evolution in the Bgc- 
M200 scaling relation, or is simply a richness-selected Ed- 
dington bias^*' is impossible to determine using only two 
clusters. Both Gilbank et al. (2007) and Andreon et 
al. (2007) found that for a small sample of clusters at 
z ^ 1 the cluster richnesses were still consistent with 
their velocity dispersions based on relations calibrated 
at lower redshift, although both parameters have large 
uncertainties in their measurements. More clusters with 
well-determined and B^^ will be needed to test if the 
cluster scaling relations at z > 1 are similar to those at 
lower redshift. 

7. SUMMARY 

We have presented a brief summary of observations for 
the northern component of the SpARCS survey. Using 
Cemini N&S spectroscopy we confirmed two rich cluster 
candidates at z ~ 1.2 selected from early survey data. 
We find that the photometric redshifts from the color 
of the cluster red-sequence agree extremely well with 
the spectroscopic redshifts. Both clusters have a smaller 
M200 than would be expected from their richness if we 
use the B^c - M200 scaling relation calibrated at z ^ 0.3. 
Whether this represents a true evolution in the cluster 
scaling relations at z > 1.2 or is simply a selection bias 
will require well-determined M200 for a larger sample of 
clusters. 

Overall, the confirmation of both SpARCS 
J163435-h402151 and SpARCS J163852-h403843 as 
bona fide massive clusters at z > 1 provide strong 
evidence that the red-sequence technique is an effective 
and efficient method for detecting clusters at z > 1 
(see also Wilson et al. 2008 who present a confirmed 
z = 1.34 cluster from the southern component of the 
SpARCS survey). The complete SpARCS catalogue 
contains hundreds of cluster candidates at z > 1 and 
promises to be one of the premier data sets for the study 
of cluster galaxy evolution at z > 1. 



Based on observations obtained with 
MegaPrime/MegaCam, a joint project of CFHT 

and CEA/DAPNIA, at the Canada- France-Hawau 
Telescope (CFHT) which is operated by the National 
Research Council (NRC) of Canada, the Institut Na- 
tional des Sciences de I'Univers of the Centre National 
de la Recherche Scientifique (CNRS) of Prance, and the 
University of Hawaii. This work is based in part on data 
products produced at TERAPIX and the Canadian 
Astronomy Data Centre as part of the Canada-France- 
Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini 

We followed up two of the richest clusters in our early dataset. 
The cluster mass fuuctiou is steep at high redshift and low mass 
systems greatly outnumber high mass systems. Due to scatter 
in the mass-richness relation lower mass systems with abnormally 
high richnesses may be more common than truly massive systems. 
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Fig. 6. — Left Panel: Galaxy velocities relative to the cluster mean velocity as a function of radius for SpARCS J163435+402151. Right 
Panel: Same as left panel but for SpARCS J163852+403843. Galaxies marked with an "x" are more likely to be near-field objects than 
members of the cluster and are not used in the computation of the velocity dispersion. 
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Fig. 7. — Left Panel: z' - 3.6^m vs. 3.6^m color magnitude diagram for galaxies at R < 550 kpc in the field of the cluster SpARCS 
J163435-I-402151. Spectroscopically confirmed cluster members and foreground/background galaxies are plotted as red and blue diamonds, 
respectively. The dotted line is the best fit line to the confirmed cluster members with the slope fixed at zero. Right Panel: Same as left 
panel but for SpARCS J163852+403843. 



partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Austrahan Research 
Council (Australia), Ministerio da Ciencia e Tecnologia 



(Brazil) and SECYT (Argentina) 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. 
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